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ABSTRACT 

Elliptical galaxies are systems where dark matter is usually less necessary to explain 
observed dynamics than in the case of spiral galaxies, however there are some instances 
where Newtonian gravity and the observable mass are insufficient to explain their ob- 
served structure and kinematics. Such is the case of NGC 4649, a massive elliptical 
galaxy in the Virgo cluster for which recent studies report a high fraction of dark mat- 
ter, 0.78 at AR e . However this galaxy has been studied within the MOND hypothesis, 
where a good agreement with the observed values of velocity dispersion is found. In a 
similar way, we have constructed a self-consistent gravitational equilibrium dynamical 
model for this galaxy assuming a modified gravity force law, which is equivalent to 
MOND for a < ao and recovers the Newtonian values for a > do- The modified gravity 
regime will be characterised by centrifugal equilibrium or dispersion velocities which 
become independent of distance, and which scale with the fourth root of the total 
baryonic mass, V A cx (MGao). We find that the recent detailed observations of the 
surface brightness profile and the velocity dispersion profile for this galaxy are con- 
sistent with the phenomenology expected in MONDian theories of modified gravity, 
without the need of invoking the presence of any hypothetical dark matter. 
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1 INTRODUCTION 

The flat rotation curves of spiral galaxies have been well 
known for decades, e.g., Bosnia (1981); Rubin et al. (1982), 
assuming the validity of the Newtonian law of gravity 
throughout, this results in a discrepancy between the dy- 
namical mass and luminous mass of spiral galaxies. Other 
types of galaxies exhibit mass discrepancies as well, the most 
remarkable is the case of dwarf spheroidal galaxies for which 
the velocity dispersion is well known and the mass inferred 
from their internal dynamics greatly exceeds the mass visible 
by factors reaching into the thousands e.g., Simon & Geha 
(2007). On the other hand, bright giant elliptical galaxies 
exhibit small mass discrepancies e.g., Romanowsky et al. 
(2003), usually these facts are commonly interpreted as the 
manifestation of the dark matter halo in which the different 
galaxies are immersed. This halo dominates the dynamic of 
extended galaxies with low surface brightness, and its pres- 
ence is less significant in more centrally concentrated galax- 
ies with larger surface brightness values, e.g., McGaugh & 
de Blok (1998). 

Alternatively, the discrepancy between observable mass 
in galaxies and their dynamics can be interpreted as a direct 
evidence for the failure of the current Newtonian and general 



relativistic theories of gravity, rather than the existence of a 
dark matter component. As examples, in spiral galaxies the 
flat rotation curve has been successful interpreted assum- 
ing a modification of Newtonian dynamics (MOND), e.g., 
Sanders & McGaugh (2002), or equivalently an extended 
Newtonian force law, e.g., Mendoza et al. (2011), the pro- 
jected surface density profiles and observational parameters 
of the local dSph galaxies are in agreement with a descrip- 
tion in MOND, e.g. Hernandez et al. (2010); McGaugh & 
Wolf (2010); Kroupa et al. (2010). Similarly, modified grav- 
ity approaches have been successful in explaining the general 
shape of the observed rotation curves of many dwarf and low 
surface brightness galaxies, e.g., Swaters et al. (2010), and 
many other astronomical systems have been interpreted in 
this context, such as globular clusters (GCs), e.g., Sollima 
& Nipoti (2010); Haghi et al. (2011); Hernandez & Jimenez 
(2012), the relative velocity of wide binaries in the solar 
neighbourhood, e.g., Hernandez et al. (2012), the infall ve- 
locity of the two component of the Bullet cluster, e.g., Moffat 
& Toth (2010) and gravitational lensing of elliptical galaxies, 
e.g., Mendoza et al. (2012). 

One generic prediction of modified gravity theories de- 
signed not to require the hypothesis of dark matter, is that 
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whenever accelerations fall below ao ~ 1.2 x 10 10 m/s 2 , a 
transition should occur away from Newtonian gravity e.g. 
MOND in Milgrom (1983), TeVeS in Bekenstein (2004), 
QUMOND, BIMOND in Zhao & Famaey (2010), the ex- 
tended Newtonian force law in Mendoza et al. (2011) or its 
relativistic version in Bernal et al. (2011). Regardless of the 
details, the modified regime will be characterized by equi- 
librium velocities for test particles orbiting within spheri- 
cal mass distributions which become constant with distance, 
and which scale with mass as V 4 tx M. 

To test the predictions of modified gravity schemes in el- 
liptical galaxies one needs to know the kinematics and mass 
distribution of these systems. Elliptical galaxies are hot sys- 
tems with little or no cold gas where determining their kine- 
matic is more difficult than for spiral galaxies. Their or- 
bital structure is thought to be the result of their evolution 
through a complex formation process. For determining their 
dynamics several tracers have been used, such as X-ray gas, 
GCs and planetary nebulae (PNe). Massive elliptical galax- 
ies are usually surrounded by hot and low density gas evi- 
dent through its X-ray glow. A tool to determine the mass 
distribution in these systems is to model their X-ray spectra 
to obtain density and temperature profiles for the gas, if we 
assume hydrostatic equilibrium for this component, we can 
then obtain the mass distribution and then an estimate of 
the brightness profile to be compared to the photometrically 
observed one, e.g., Das et al. (2010). 

The problem of modeling the luminosity profile and 
kinematics of elliptical galaxies has been treated before, e.g., 
Kormendy et al. (2009) modelled the brightness profiles of 
all known elliptical galaxies in the Virgo cluster using Sersic 
profiles logl cx r ,1//n , they develop a mechanism to calcu- 
late a realistic error in the Sercic parameters and identify 
departures from these profiles that are diagnostic of galaxy 
formation processes. Teodorescu et al. (2011) studied the 
kinematics of PNe in the Virgo giant elliptical galaxy NGC 
4649 (M 60), assuming Newtonian gravity, they conclude 
that the kinematics of this object are consistent with the 
presence of a dark matter halo around M60, this halo is al- 
most one-half of the total mass of the galaxy within 3R e ■ In 
De Bruyne et al. (2001) three-integral axisymmetric mod- 
els for NGC 4649 and NGC 7097 are considered, concluding 
that the kinematic data of NGC 4649 are consistent with 
a dynamical model with a moderate amount of dark mat- 
ter, Das et al. (2010) create a dynamical model of NGC 
4649 using the NMAGIC code and kinematic constraints to 
infer a dark matter mass fraction in NGC 4649 of ~ 0.78 
at 4Re. In Samurovic & Cirkovic (2008a) and Samurovic & 
Cirkovic (2008b) the GC dynamics of NGC 4649 are used 
as mass tracers and the Jeans equation is solved for New- 
tonian and MOND models, finding that both are consistent 
with the values of the observed velocity profile of the galaxy, 
although considering the light distribution only to first ap- 
proximation, as a radial power law for the surface brightness 
profile. 

Following this ideas, in this paper we construct fully 
self-consistent dynamical models for the massive elliptical 
galaxy NGC 4649 using the modified Newtonian force law 
formulation of Mendoza et al. (2011), to explore the consis- 
tency of that approach. The free parameters of the galactic 
model are calibrated to match the details of the observed 
surface brightness and projected velocity dispersion profiles. 



We find that fully self-consistent equilibrium models can 
be constructed to match all observational constraints avail- 
able on NGC 4649, including the outward flattening of the 
projected velocity dispersion profiles, without invoking any 
dark matter. To within uncertainties, the same happens un- 
der Newtonian gravity, but requiring substantial amounts of 
the hypothetical and as yet undetected dark matter. 

In section (2) we derive the model through which equi- 
librium profiles for spherically symmetric stellar populations 
are derived, under the modified Newtonian force law of Men- 
doza et al. (2011). In section (3) we show that such models 
can be easily obtained to satisfy all observed parameters for 
NGC 4649, all showing clearly a flattening of their projected 
velocity dispersion profiles at large radius. Our conclusions 
appear in section (5). 



2 NON-ISOTHERMAL GRAVITATIONAL 
EQUILIBRIUM MODELS 

In the same way as done for globular clusters e.g. Sollima 
& Nipoti (2010), in Hernandez & Jimenez (2012) we have 
modeled the massive elliptical galaxy NGC 4649 as a popula- 
tion of self gravitating stars in a spherically symmetric equi- 
librium configuration, under a modified Newtonian gravita- 
tional force law. For a test particle at a distance r from the 
centre of a spherically symmetric mass distributions M(r): 

In the above x = l M /r, where l M = (GM(r)/a ) 1/2 . We 
see that when a >> ao, x >> 1 and f(x) — > aox 2 , one 
recovers Newton's force law, while when a « ao, x « 1 
and f(x) — s- aox, where an equivalent MOND force law / = 
(M(r)Gao) 1 / 2 /r is obtained. The index n > 1 mediates the 
abruptness of the transition between the two regimes. 

This was shown in Mendoza et al. (2011) to yield gen- 
eralised isothermal gravitational equilibrium configurations 
with characteristic finite radii, masses and velocity disper- 
sion values, r g , M, a, which smoothly evolve from the clas- 
sical virial equilibrium of M = a 2 r g /G to the observed 
tilt in the fundamental plane of elliptical galaxies, to the 
a 4 = (MGao) scaling of the galactic Tully-Fisher relation, 
as one goes from x >> 1 to x ~ 1 to x « 1. Also, con- 
sistency with solar system observations was found there to 
constrain the transition to be fairly abrupt, requiring n > 4, 
and in Mendoza et al. (2011) it is also proven that a sufficient 
condition for Newton's theorems for spherically symmetric 
mass distributions to hold, for any modified force law, is 
only that / can be written as a function exclusively of the 
variable x. 

Recently, in Hernandez & Jimenez (2012) we showed 
that constructing similar successful models for globular clus- 
ters requires n > 10 to obtain best fits to brightness and 
velocity dispersion profiles, here we use n = 10 to construct 
equilibrium models for the elliptical galaxy NGC 4649. No- 
tice also that given the form of the proposed force law, any 
other larger value of n results in only marginal differences 
from n = 10. In addition, Bernal et al. (2011) showed this 
modified force law model to be the low velocity limit of a 
formal generalization to GR of the f(R) type, providing a 
theoretical basis for the model used. 
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Taking the derivative of the kinematic pressure, the 
equation of hydrostatic equilibrium for a polytropic equa- 
tion of state P = Kfp is: 



d(Kp~> 
dr 



- P v<t>. 



(2) 



In going to locally isotropic Maxwellian conditions, 7=1 
and K = er(r) 2 and we get: 



dr 



dr 



-V0. 
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since p = (47rr 2 ) 1 dA J^ the preceding equation can be 



written as: 
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where a(r) is the isotropic Maxwellian velocity dispersion for 
the population of stars, which is allowed to vary with radius, 
as observed in NGC 4649, e.g., Pinkney et al. (2003). The 
above is a generalization of the treatment presented in Her- 
nandez et al. (2010), which we used in the modeling of dSph 
galaxies, which are characterized by flat velocity dispersion 
profiles, obtaining mass models consistent with observed ve- 
locity dispersion, half mass radii and total masses, in the 
absence of dark matter. Another example of a similar treat- 
ment is found in Hernandez & Jimenez (2012), where we 
modeled eight GCs in our Galaxy proposing a function a(r) 
for modeling the velocity dispersion profile of GCs assum- 
ing the same modified law of gravity used here, accurately 
reproducing all observational constraints. 

As an illustrative example we can take the limit f(x) = 
aox 2 , and recover —GM(r)/r 2 for the right hand side of 
equation (3), the Newtonian expression appearing for a >> 
ao, or equivalently x « Im. If one then imposes isothermal 
conditions er(r) = a and looks for a power law solution for 
M(r) = M (r/r ) m , we get: 



GMo 



(5) 



and hence m = 1, the standard isothermal halo, M(r) — 
2a 2 r/G, having a constant centrifugal equilibrium velocity 
v 2 = 2a 2 and infinite extent. At the other limit, a << ao, 
x >> Im, equation (3) yields: 



[GM a ] 1/2 



(- 

\r 



1/2 



(6) 



In this limit m — 0, we obtain M(r) = Mo and v 2 = 3a 2 = 
(GMoao) 1 ^ 2 , the expected Tully-Fisher scaling of the cir- 
cular equilibrium velocity with the fourth root of the mass, 
with rotation velocities which remain flat even after the mass 
distribution has converged, rigorously isothermal halos are 
naturally limited in extent. It is interesting to note that in 
this limit the scaling between the circular rotation velocity 
and the velocity dispersion is only slightly modified as com- 
pared to the Newtonian case, the proportionality constant 
changes from 2 to 3, for the squares of the velocities. In astro- 
physical units, this low acceleration limit for an isothermal 
halo in gravitational equilibrium yields: 
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Figure 1. Volumetric density profile for the best fit model for 
NGC 4649. The vertical line gives the point where x = 1. Note 
that the asymptotic p(r) profile at large radii is steeper than the 
r~ 2 of the Newtonian case, resulting in a finite total radius and 
total mass for the configuration. 
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km/s, 



(7) 

(8) 



for the velocity dispersion and the centrifugal equilibrium 
velocities of a halo of total mass Mo, respectively. 

Taking initial conditions M(r) — > 0, dM(r)/dr — > 
Anr 2 po, when r — > 0, a constant central density po, we can 
solve the full second-order differential equation (4) for M(r) 
through a numerical finite differences scheme, once a model 
for a(r) is adopted. 

Solving equation (4) yields the volumetric profiles for 
the density and mass, p(r), M(r). p(r) is then projected 
along one dimension to obtain a projected surface density 
mass profile, T,(R). In all that follows we shall use r for ra- 
dial distances in 3D, and R for a projected radial coordinate 
on the plane of the sky. E(ii) can be compared to observed 
surface brightness profiles once a mass-to-light ratio is as- 
sumed. 

Finally, a projected velocity dispersion profile cr p (R) is 
constructed through a volume density weighted projection 
of the volumetric velocity dispersion profile of a model, cr(r). 
Therefore, a proposed model for cr(r) does not directly give 
a projected cr P (R) profile, which the model only yields after 
having solved for the detailed density structure, and the 
subsequent mass weighted projection of the proposed cr(r). 

Observations provide only the density weighted pro- 
jected profiles of a(r), integrated along the line of sight, 
<jp(R), not the volumetric u(r) profiles which we require in 
equation (4). We shall therefore adopt a parametric form 
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Figure 2. Volumetric mass profile for the best fit model for NGC 
4649. The vertical line gives the point where x = 1. Note that the 
total mass converges for a finite radius. 

for the volumetric velocity dispersion profile, and adjust the 
parameters to obtain a match to the observed galactic prop- 
erties. For this we use: 

a(r) = CTiexp ) + °oo- (9) 

In the above <Too is given directly by the observations 
as the asymptotic value of the measured projected velocity 
dispersion profile for the galaxy, this value is ~ 200fcms~ 1 , 
Das et al. (2011), as at larger R radii, projection effects 
tend to zero and o-(R) — > o{r). This leaves us with three 
model parameters to determine: po, ci and r CT , which we 
fit to match the observed projected velocity dispersion pro- 
file, as well as a central projected velocity dispersion value 
for NGC 4649 of » AOOKms' 1 Pinkney et al. (2003), and 
comparing the resulting model £(i?) profile to the observed 
surface brightness profile of NGC 4649, taking M/L = 8.0, 
following Bridges et al. (2006). 



3 MODELLING ELLIPTICAL GALAXY NGC 
4649 

We begin this section presenting in figures (1) and (2) our 
best fit NGC 4649 model. In figure (1) we show the volumet- 
ric density profile, which is qualitatively similar to a cored 
isothermal profile, with the difference that the asymptotic 
p(r) profile at large radii is steeper than the r~ 2 of the New- 
tonian case. This results in finite total masses and finite half- 
mass radii even for the asymptotically flat cr(r) volumetric 
profiles we adopt, in contrast with the situation in classical 
gravity, where infinitely extended mass profiles would ap- 
pear. This finite profiles are also what appears under the 
type of modified gravity laws we are treating, even for rig- 



Figure 3. Comparison of the resulting model projected surface 
brightness for NGC 4649 (thick continuous line) to the corre- 
sponding observed quantities and the best fit S'ersic function 
with the central (dashed line), maximum and minimum values 
consistent with the reported confidence intervals for the S'ersic 
parameters, flanking curves. The squares correspond to data in 
(Kormendy et al. 2009), and the points with error bars for PNe 
data, (Tcodorcscu et al. 2011) Following Bridges et al. (2006) 
we take M/Ly = 8, which coincides with the stellar population 
studies on Shen & Gebhardt (2010). The vertical line gives the 
point where x = 1 and the modified force law used shifts from the 
classical Newtonian form to the MONDian character of eq. (6) 



orously isothermal a(r) = ero equilibrium configurations, as 
already pointed out in Hernandez et al. (2010) and Mendoza 
et al. (2011) and as shown in the developments following eq. 
(6). The vertical line shows the radius where the condition 
x = 1 is crossed, i.e., the point beyond which a < ao, and 
hence the departure of the force law used from the New- 
tonian value towards the MONDian regime, consequently 
the point where the density profile steepens away from the 
1/r 2 of Newtonian gravity, towards the convergence of the 
MONDian regime. 

Figure (2) shows the corresponding volumetric radial 
mass profile for the same model, indicating again with a ver- 
tical line the threshold where x = 1. This clearly shows the 
departure from the linear growth of the total mass within the 
approximately isothermal Newtonian region, towards con- 
vergence to a finite total mass in the MONDian a < ao 
regime. We see that 40% of the total model mass lies be- 
yond x = 1. This fraction is distributed over a much larger 
area having much smaller projected mass surface densities 
than the central regions interior to the x — 1 threshold, the 
ones which are much more easily observed, and over which 
Newtonian gravity accurately holds. 

These two figures illustrate the physical consistency of 
the model, a positive isotropic Maxwellian volumetric veloc- 
ity dispersion is assumed, integration of eq.(4) then yields 
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Figure 4. Comparison of the resulting model projected velocity 
dispersion profile for NGC 4649 (continuous line) to the corre- 
sponding observed quantities, stars for long slit observations in 
(Pinkncy et al. 2003) and the points with error bars for PNe data. 
The vertical line gives the point where a = ao- 

a volumetric density profile essentially consistent with what 
would appear under Newtonian dynamics in the region in- 
terior to x = 1, where the force law converges precisely to 
the standard expression. At large radii however, the density 
profile increasingly steepens and naturally reaches p(r) = 
at a well denned total radius, as is evident from the conver- 
gence seen in figure (2). Thus, the distribution function is 
necessarily positive throughout the modeled structure, and 
goes to zero at a well denned outer radius. 

The best fit model constructed as described above was 
optimised to match the inferred total mass, the observed 
projected surface brightness profile under the assumption of 
constant M/Lv, and the projected velocity dispersion mea- 
surements. We take the surface brightness profile observa- 
tions from photometry in Kormendy et al. (2009) and the 
scaled number density calculated from PNe in Teodorescu 
et al. (2011) as presented in Das et al. (2011), and projected 
velocity dispersion measurements from Pinkney et al. (2003) 
adapted in De Bruyne et al. (2001). Following Bridges et al. 
(2006) we take M/Lv = 8, which coincides with the cen- 
tral value reported by Shen & Gebhardt (2010) determined 
directly from stellar population synthesis models. 

The final resulting surface brightness profile and its cor- 
responding volume density weighted projected velocity dis- 
persion profile <J P (R) for galaxy NGC 4649 are presented in 
figures 3 and 4, the parameters for the best fit in eq.(9) are 
o~oo = 175km/ s, (Ti = 230fcm/s and r a = 8.5fcpc, the central 
density that we have assumed is po = 1-5 x 10 2 M Q /pc 3 . 
In the figures the vertical line indicates the point where 
a — ao, we can see that the model accurately fits the ob- 
served projected surface brightness profiles of the galaxy and 
the central value of the measured projected velocity disper- 



sion o p (R = 0) = 400fcm/s, as well as the observed profiles 
for <T P (R). 

We see that both, the surface brightness and velocity 
dispersion profiles are in agreement with photometry and 
PNe observations, to within reported uncertainties. Under 
Newtonian gravity this requires to invoke a large amount 
of dark matter to reproduce the structure and kinematics 
of this galaxy, while under the modified force law scheme 
treated here, it is a natural consequence of the change in 
gravitational regime, leading to finite matter distributions 
when a < ao and a(r) tends to a constant. 

The total mass we obtain from the dynamical model- 
ing described is consistent with what results from the in- 
tegration of the surface brightness profile assuming a con- 
stant mass to light ratio M/Lv = 8. The mass obtained 
for the dynamical model to NGC 4649 is 1.09 x 10 12 M Q , 
and the luminosity that results from integrating the V 
band light profile is 1.36 x W 11 Lq, implying a mass of 
1.08 x lO 12 Af from using M/Lv = 8 as determined di- 
rectly from the stellar population studies of Shen & Geb- 
hardt (2010). Also, notice that the asymptotic value of the 
velocity dispersion is 195.0±30.36 km/s in consistency with 
the a = 0.2(M/M Q ) 1/4 = 203.9fcm/s prediction of eq. (7) 
for 1.08 x 10 12 M Q . 



4 CONCLUSIONS 

We show that for the massive elliptical galaxy NGC 4649, 
fully self consistent spherically symmetric equilibrium mod- 
els can be constructed using a modified Newtonian force 
law which smoothly transits from the Newtonian value when 
a > ao, towards the MONDian phenomenology in the a < ao 
regime, which naturally satisfy all observational constraints 
available for the central, asymptotic and radial profiles for 
projected surface brightness and velocity dispersion mea- 
surements. All this, without the need of invoking any hypo- 
thetical and as yet undetected dark matter component. 

The resulting models are typically characterized by a 
Newtonian inner region, which smoothly transits to a mod- 
ified gravity outer region on crossing the a = ao thresh- 
old. The resulting internal velocity dispersion profiles corre- 
spondingly transit from an inner radially decaying region to 
an outer flat velocity dispersion one. 

We show that the asymptotic value of the measured 
velocity dispersion profile, a p (R — > oo), and the total mass 
for this systems, M, are consistent with the generic modified 
gravity prediction for a p (R -> oo) = O.2(M/M ) 1/4 . 
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